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Abstract 

\o : 

ON ■ 

We consider the Hamiltonian dynamics and thermodynamics of spherically 
Q H . symmetric Einstein-Maxwell spacetimes with a negative cosmological con- 

$h . stant. We impose boundary conditions that enforce every classical solution to 

be an exterior region of a Reissner-Nordstrom-anti-de Sitter black hole with 
' a nondegenerate Killing horizon, with the spacelike hypersurfaces extending 

^ ■ from the horizon bifurcation two-sphere to the asymptotically anti-de Sit- 

ter infinity. The constraints are simplified by a canonical transformation, 
which generalizes that given by Kuchaf in the spherically symmetric vacuum 
Einstein theory, and the theory is reduced to its true dynamical degrees of 
freedom. After quantization, the grand partition function of a thermody- 
namical grand canonical ensemble is obtained by analytically continuing the 
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Lorentzian time evolution operator to imaginary time and taking the trace. 
A similar analysis under slightly modified boundary conditions leads to the 
partition function of a thermodynamical canonical ensemble. The thermo- 
dynamics in each ensemble is analyzed, and the conditions that the (grand) 
partition function be dominated by a classical Euclidean black hole solution 
are found. When these conditions are satisfied, we recover in particular the 
Bekenstein-Hawking entropy. The limit of a vanishing cosmological constant 
is briefly discussed. 

Pacs: 04.60.Ds, 04.60.Kz, 04.70.Dy, 04.20.Fy 
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I. INTRODUCTION 



Hawking's celebrated result of black hole radiation |TJ and related developments 
made it possible to consider thermodynamical equilibrium systems involving black holes in 
the manner first anticipated by Bekenstein |5|||. At a semiclassical, "phenomenological," 
level, a black hole thermodynamical equilibrium system can be introduced by simply im- 
mersing a radiating black hole in a heat bath such that the outgoing Hawking radiation 
balances the radiation that falls in from the bath ]|7|-|TT|j. At a deeper level, one aspires to 
construct a full thermodynamical equilibrium ensemble by starting from a quantum the- 
ory of gravity for black hole type geometries j|T2"|-|T7|] . For reviews, see for example Refs. 



16.18H20 



At the semiclassical level, the thermodynamical equilibrium configurations involving 
black holes tend to be unstable against thermal fluctuations The classic example 

is a Schwarzschild black hole in equilibrium with an asymptotically flat heat bath, in the 
approximation where the back-reaction of the radiation on the geometry is neglected: the 
heat capacity in this instance is — (871T 2 ) , where T is the temperature measured at the 
infinity, and the fact that this heat capacity is negative indicates thermodynamical insta- 
bility. While such instabilities are not unexpected in self-gravitating systems, they do pose 
an obstacle to constructing thermodynamical equilibrium ensembles from quantum gravity. 
This is because the existence of a thermodynamical ensemble implies the positivity of cer- 
tain response functions associated with that ensemble PTf . For example, in the canonical 
ensemble the heat capacity is necessarily positive; consequently, a canonical ensemble of the 
usual kind does not appear to exist for Schwarzschild black holes in asymptotically flat space 



22]. 

To construct a thermodynamical ensemble appropriate for black hole geometries from 
a quantum theory of gravity, one thus needs to choose the boundary conditions for the 
ensemble in a judicious manner, motivated by the stability of the corresponding semiclassical 
equilibrium situations. One possibility is to replace an asymptotic infinity by a finite "box" 



at which the local temperature is then fixed [14,23-34]. The possibility on which we shall 



concentrate in this paper is to include a negative cosmological constant Jl 

A negative cosmological constant makes classical black hole solutions asymptotically anti- 
de Sitter. We shall consider spherically symmetric spacetimes, and as the only matter field 
we include the spherically symmetric Maxwell field. All the relevant classical solutions then 
belong to the Reissner-Nordstrom-anti-de Sitter (RNAdS) family p0|-43]. The temperature 
of the Hawking radiation is redshifted to zero at the asymptotically anti-de Sitter infinity, 
but from the rate at which the local Hawking temperature approaches zero one can extract a 
"renormalized" Hawking temperature, and this renormalized Hawking temperature can then 



be taken as one fixed quantity in the thermodynamical ensembles |lq , |35| , |36| , |38| . We shall 



consider both the canonical ensemble, in which the electric charge is fixed, and the grand 
canonical ensemble, in which the electric potential difference between the event horizon and 
the infinity is fixed. 

To quantize the theory and to build the equilibrium ensembles, we shall adapt the method 
introduced in Ref. |$2| in the context of spherically symmetric vacuum geometries in the 
presence of a finite boundary. We shall first set up a classical Lorentzian Hamiltonian 
theory in which, on the classical solutions, the right end of the spacelike hypersurfaces is 
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at the asymptotically anti-de Sitter infinity in an exterior region of a black hole spacetime, 
and the left end of the hypersurfaces is at the bifurcation two-sphere of a nondegenerate 
Killing horizon. We then canonically quantize this theory, and obtain the thermodynamical 
(grand) partition function by suitably continuing the Schrodinger picture time evolution 
operator to imaginary time and taking the trace. A crucial input is how to handle the 
analytic continuation at the bifurcation two-sphere. As in Ref. | 32|j , we shall see that a 
continuation motivated by smoothness of Euclidean black hole geometries yields a (grand) 
partition function that is in agreement with that obtained via path integral methods. 

To implement the method used in Ref. [ 32" | , one must be able to canonically quantize the 
Lorentzian theory in some practical fashion. In Ref. [|32] this was achieved by using canon- 
ical variables that were first introduced by Kuchaf under asymptotically flat, Kruskal-like 
boundary conditions MM. In these variables the constraints of the vacuum theory become 



exceedingly simple, and the classical Hamiltonian theory can be explicitly reduced into an 
unconstrained Hamiltonian theory with just one canonical pair of degrees of freedom. We 
shall show that an analogous set of canonical variables exists for our system, and the classi- 
cal Hamiltonian theory can again be explicitly reduced into an unconstrained Hamiltonian 
theory. Under boundary conditions tailored to the grand canonical ensemble, the reduced 
Lorentzian Hamiltonian theory has two canonical pairs of degrees of freedom;[| under bound- 
ary conditions tailored to the canonical ensemble, the reduced Lorentzian Hamiltonian the- 
ory has just one pair of canonical degrees of freedom. Using these variables, it will be 
possible to construct a quantum theory and a (grand) partition function in close analogy 



with Ref. 32 



It will turn out that both the canonical ensemble and the grand canonical ensemble for our 
system are well defined. In particular, the appropriate thermodynamical response functions 
are positive. We shall also be able to give the conditions under which the (grand) partition 
function is dominated by a classical Euclidean solution. The grand canonical ensemble 
exhibits a transition from a region where a classical solution dominates to a region where 
no classical solution dominates, in close analogy with what happens with the spherically 



symmetric boxed vacuum canonical ensemble [24,25 . As in Refs. [24. 251, one may see this 



as evidence for a phase transition between a black hole sector and a topologically different 
"hot anti-de Sitter space" sector. In the canonical ensemble we find evidence for this kind 
of a phase transition only in the special case when the charge vanishes. When the charge 
is nonvanishing, there occurs a different kind of phase transition in which the dominating 
contribution to the partition function shifts from one classical solution to another as the 
boundary data changes. 

The rest of the paper is as follows. In Sec. |TJ we set up a classical Hamiltonian theory un- 
der boundary conditions tailored to the grand canonical thermodynamical ensemble, paying 
special attention to the falloff conditions at the asymptotically anti-de Sitter infinity |f46| , [47 



In particular, we choose to fix the values of the electric potential at the infinity and at the 



1 The conclusion of two canonical pairs of degrees of freedom for the spherically symmetric Einstein- 
Maxwell system with a cosmological constant was previously reached, under a different set of 
boundary conditions, in Ref. [45|. 
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horizon, in a manner that will be made precise in terms of the fall-off conditions. In Sec. 
1 1 1 we perform the canonical transformation, and in Sec. IV the constraints are eliminated 
and the theory is reduced to its true dynamical degrees of freedom. In Sec. |V] we quantize 
the theory and obtain the grand partition function of the thermodynamical grand canonical 
ensemble. The thermodynamics in the grand canonical ensemble is analyzed in Sec. [VT|. 
Sec. [VIl] outlines the corresponding classical, quantum mechanical, and thermodynamical 
analyses under boundary conditions that fix the charge instead of the electric potential, and 
thus lead to the thermodynamical canonical ensemble. 

The results are summarized and discussed in Sec. |V111| . Some facts about the RNAdS 
solutions are collected in Appendix Finally, Appendix [B] outlines the classical Hamil- 
tonian analysis and the quantization of the reduced Hamiltonian theory in the case where 
the cosmological constant vanishes and the asymptotically anti-de Sitter falloff conditions 
are replaced by asymptotically flat falloff conditions. With asymptotic flatness, neither the 
partition function nor the grand partition function turns out to be well-defined, and we 
recover neither a canonical ensemble nor a grand canonical ensemble. 



II. CANONICAL FORMULATION IN THE METRIC VARIABLES 

In this section we present a Hamiltonian formulation of spherically symmetric electrovac- 
uum spacetimes with a negative cosmological constant, with boundary conditions appropri- 
ate for the exterior of a RNAdS black hole with a nondegenerate horizon. Some relevant 
properties of the RNAdS metric are reviewed in Appendix [A]. 

We consider the general spherically symmetric ADM metric 

ds 2 = -N 2 dt 2 + A 2 (dr + N r dt) 2 + R 2 dQ 2 , (2.1) 

where dfl 2 is the metric on the unit two-sphere, and N, N r , A and R are functions of t and 
r only. The electromagnetic potential is taken to be described by the spherically symmetric 
one-form 

A = Tdr + <&dt , (2.2) 

where V and $ are functions of t and r only. The fact that this one-form is globally 
defined makes the electromagnetic bundle trivial, and will preclude the black hole from 
having a magnetic charge. The coordinate r takes the semi-infinite range [0, oo). Unless 
otherwise stated, we shall assume both the spatial metric and the spacetime metric to be 
nondegenerate. In particular, A, R, and N are taken to be positive. We shall work in natural 
units, h — c — G — 1. 

The action of the Einstein-Maxwell theory with a negative cosmological constant is 

S = — fd 4 x J-Wg ( {i) R + U~ 2 - F^Fj) + (boundary terms) , (2.3) 

167T J ^ ' 

where is the determinant of the four-dimensional metric, ^R is the four- dimensional 
Ricci scalar, and F^ v = d fl A u — d u A^ is the electromagnetic field tensor. The cosmological 
constant has been written as — 3£~ 2 , where i > 0. Inserting the spherically symmetric fields 
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( j2.ip and ( p.2|) and integrating over the two-sphere we obtain, up to boundary terms, the 
action 



SxiA,R,T;N,N r ,<S>] 



dt 



dr 



-N' 1 (R(-A + {AN r )'){-R + R'N r ) + U(-i? + R'N r 



+^N- l A- 1 R 2 (t - $') 

+N (A- 2 RR'A' - A^RR" - \A- l R' 2 + \A + \r 2 AR 2 ) 



(2.4) 



The equations of motion derived from local variations of ( |2.4j ) are the full Einstein-Maxwell 
equations for the spherically symmetric fields ( |2.1| ) and ( |2.2|) . A generalized Birkhoff's 
theorem can be proven using the same techniques as in the case of a vanishing cosmological 
constant ||48|| : every classical solution is locally either a member of the extended RNAdS 
family (see Appendix [A]), or a spacetime that generalizes the Bertotti- Robinson solution to 
accommodate a negative cosmological constant [p|- ^3| , |48[] . We shall address the boundary 
conditions and boundary terms that are needed to make the variational principle globally 
well-defined after passing to the Hamiltonian formulation. 

The momenta conjugate to the configuration variables A, R, and T are 



P A = -N^RiR - R'N r ) 
P R = -AT- 1 (A(R - R'N r 

P r = N^A^R 2 ^ - $') 



R(A - (AN r )'] 



(2.5a) 
(2.5b) 
(2.5c) 



A Legendre transformation gives the Hamiltonian action 



Sx[A,R,r,P A ,P R ,P r ;N,N r ,$) 



dt 



dr ( P A A + P R R + P r T -NH- N r H r - $G 



(2.6) 



where the super-Hamiltonian constraint H, the radial supermomentum constraint H r , and 
the Gauss law constraint G are given by 



H 



G 



R-'PrPa 



\R- 2 A{P 2 + P 2 ] 



+A- 1 RR" - A- 2 RR'A' + ^A^R' 2 



P R R ! - AP A 



|A - f t- 2 AR 2 



TP* , 



-P' 



(2.7a) 
(2.7b) 
(2.7c) 



We have written the electric potential $ in terms of the quantity 

$ := $ - N r T , 



(2.8) 



which now acts as the Lagrange multiplier associated with the Gauss constraint in ( ^.6[ ). 
It would be possible to proceed retaining $ as the Lagrange multiplier, and the super- 
momentum constraint would then be the same as without the electromagnetic field (see, 
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for example, Ref. |49|]). However, using <3> has the technical advantage that the supermo- 
mentum constraint ( |2.7b|) generates spatial diffeomorphisms in both the gravitational and 
electromagnetic variables. This fact will prove useful in Sec. fT\\ . 

The Hamiltonian equations of motion are obtained from local variations of (|2.6|) . The 
constraint equations are 



H 

H r 
G 









and the dynamical equations of motion read 

A = N{R~ 2 AP A - R- X P R ) + (iV r A)' , 

R = -NR^Pa + N r R' , 

f = NAR~ 2 P T + (N r V)' + !>' , 

2 , - . n -2 D 2 



Pa 
Pr 



2 

N 



N r Pl 



\N [-R- 2 {Pl + P r 2 ) - (A^R'Y + 1 + 3£-' z R 

~AR~ 3 (Pl + P 2 ) - R~ 2 P A P R - (A~ L R')' + 3£~' 2 AR 
(A-'N'Ry + i^Pn)' , 
r • 



- A- 2 N'RR' + N r P' K , 



(2.9a) 
(2.9b) 
(2.9c) 



(2.10a) 
(2.10b) 
(2.10c) 
(2.10d) 

(2.10e) 
(2.10f) 



It is easy to verify that the Poisson bracket algebra of the constraints closes, and we thus 
have a first class constrained system |50 . 



We now wish to adopt boundary conditions that enforce every classical solution to be 
an exterior region of a RNAdS spacetime with a nondegenerate horizon (see Appendix |A|), 
such that the constant t hypersurfaces begin at the horizon bifurcation two-sphere at r = 
and reach the asymptotically anti-de Sitter infinity as r — >• oo. 

Consider first the left end of the hypersurfaces. At r — > 0, we adopt the conditions 



A(t,r) 


= A (t) + 0(r 2 ) , 


(2.11a) 


R(t,r) 


= R (t)+R 2 (t)r 2 + O(r A ) , 


(2.11b) 


PA{t,r) 


= 0(r 3 ) , 


(2.11c) 


P R (t,r) 


= 0(r) , 


(2.11d) 


N(t,r) 


= Aq(t)r + 0(r 3 ) , 


(2. lie) 


N r (t,r) 


= N{(t)r + 0(r 3 ) , 


(2.11f) 


T{t,r) 


= 0{r) , 


(2.11g) 


P r (t,r) 


= Qo(t) + Q 2 (t)r 2 + 0(r 4 ) , 


(2.11h) 


$(£, r) 


= $ (t) + O(r 2 ) , 


(2-Hi) 



where Ao and Rq are positive, and N\ > 0. Here 0(r n ) stands for a term whose magnitude 
at r — > is bounded by r n times a constant, and whose fc'th derivative at r — > is similarly 
bounded by r n_fc times a constant for 1 < < n. It is straightforward to verify that these 
falloff conditions are consistent with the constraints H = H r = G = 0, and that they are 
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preserved by the time evolution equations. The metric falloff conditions ( 2.11a|) - (|2.11f ), 



which are identical to those introduced in Ref. [32[ in the context of the Schwarzschild 
black hole, guarantee that the classical solutions have a nondegenerate horizon, and that 
the constant t hypersurfaces begin at r = at a horizon bifurcation two-sphere in a manner 
asymptotic to hypersurfaces of constant Killing time.Q The coordinates become thus singular 
at r — > 0, but this singularity is quite precisely controlled. In particular, on a classical 
solution the future unit normal to a constant t hypersurface defines at r — > a future 
timelike unit vector n a (t) at the bifurcation two-sphere, and the evolution of the constant t 
hypersurfaces boosts this vector according to 

n a (tx)n a {t 2 ) = - cosh QT a A^^N^t) dtj . (2.12) 

The falloff conditions ( |2.11gD -( [2.11i|) for the electromagnetic field variables are motivated 
by our thermodynamical goal, and they will be discussed further in Sec. |V|. 

Consider then the right end of the hypersurfaces. At r —>■ oo, we assume that the variables 
have asymptotic expansions in integer powers of (1/r), with the leading order behavior 



A{t,r) 


= ir- 1 - §fV~ 3 + \(t)£ 3 r- A + 0°°(r- 5 ) , 


(2.13a) 


R(t,r) 


= r + fp{t)r- 2 + O co {r- 3 ) , 


(2.13b) 


PA(t,r) 


= 0°°(r- 2 ) , 


(2.13c) 


P R (t,r) 


= 0°°(r- 4 ) , 


(2.13d) 


N(t,r) 


= A^R' (N + (t) +0°°(r~ 5 )) , 


(2.13e) 


N r {t,r) 


= 0°°(r- 2 ) , 


(2.13f) 


F(t,r) 


= 0°°(r- 2 ) , 


(2.13g) 


P r {t,r) 


= Q + (t) + 0°°(r- 1 ) , 


(2.13h) 


$(t,r) 


= $ + (t) + 0° o (r- 1 ) , 


(2.131) 



where N+(t) > 0. 0°°(r~ n ) denotes a term that falls off at infinity as r~ n , and whose deriva- 
tives with respect to r fall off accordingly as r~ n ~ k , k = 1, 2, . . .. It is again straightforward 
to verify that these falloff conditions are consistent with the constraints and that they are 
preserved by the time evolution equations. Comparison with Ref. |47j shows that the met- 



ric is asymptotically anti-de Sitter, with the constant t hypersurfaces being asymptotic to 
hypersurfaces of constant Killing time, and N+{t) gives the rate at which the Killing time 
evolves with respect to t at the infinity. Note that the lapse-function N diverges at the 
infinity for any nonzero value of N(t). For future use, we define the quantity 

M+(t) := A(t) + 3p(t) . (2.14) 



2 The text in Ref. [ 32 j contains at this point a minor inaccuracy. Equations (2.6a) and (2.6b) of Ref. 
[32] [our ( |2.11a| ) and (2.11b)] are not sufficient to ensure that the hypersurfaces end at the horizon 
bifurcation two-sphere, but for example the set (2.6a)-(2.6c) of Ref. [ 32 1 [our ( 2.11a[) -( ^.llc )1 is. 
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When the equations of motion hold, M+(t) is independent of t, and it is equal to the mass 
parameter of the RNAdS metric (|Al|). 

Taken together, the falloff conditions ( 2.11|) and ( p,13| ) achieve our aim. Every classical 



solution is an exterior region of a RNAdS spacetime with a nondegenerate event horizon, 
such that the constant t hypersurfaces begin at the horizon bifurcation two-sphere and 
reach the asymptotically anti-de Sitter infinity. In particular, the classical solutions satisfy 
i? 2 >0.Q 

It would be possible to replace (p. 13a] ) and ( |2.13b| ) by 



A(t,r) =£r- 1 -|£ 3 r- 3 + 0°°(r- 4 ) , (2.15a) 
R(t,r) = r + 0°°(r- 2 ) , (2.15b) 

and then drop the assumption that the expansion proceed in integer powers of (1/r) beyond 
the order shown, provided one makes more precise assumptions about what is meant by the 
symbol 0°°. Alternatively, it would be possible to strengthen the falloff conditions to read 

A(t,r) = fr" 1 - i£ 3 r- 3 + A(t)£ 3 r~ 4 + 0°°(r- 4 - e ) , (2.16a) 
R(t,r) = r + 0°°(r~ 2 - t ) , (2.16b) 



where < e < 1, with similar changes in the rest of (|2.13 ). This would be analogous to 



the falloff conditions adopted for the asymptotically flat Schwarzschild case in Ref. [44 



in that the value of the mass ( |2.14 ) could then be read solely from the expansion ([2.1 6a|) 



of A. One might also consider writing the theory in terms of a lapse-function that has been 
rescaled by the factor A~ 1 R': by ( |2.13e ), the falloff of the new lapse at r — > oo would then 



be independent of the canonical variables. For concreteness, we shall adhere to the theory 
as written above. 

We can now write an action principle compatible with our falloff conditions. Consider 
the total action 

S[A, R, T, P A , P R , P T ; N, N r , $] = S S [A, R, V, P A , P R , P T ; N, N r , $] 

+S d x[A,R,Q ,Q + ;N,$ ,$ + ] , (2.17) 

where the boundary action is 

S a v[A,R,Q ,Q+;N, $ ,$ + ] 

= Jdt (l^A^iA- 1 - N+M+ + $ Qo - $ + g + ) . (2.18) 



The total action ( 2.17|) is clearly well-defined under our boundary conditions. Its variation 



contains a volume term proportional to the equations of motion, boundary terms from the 



3 The falloff conditions ( [2.11 ) are compatible with either sign of i?2- The case R2 < would 
correspond to the bifurcation two-sphere of an inner horizon, which is excluded from the classical 
solutions only after the asymptotically anti-de Sitter falloff has been invoked at r — > 00. If desired, 
the requirement R2 > could of course be already added to the conditions ([2.11]) . 
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initial and final hypersurfaces proportional to SA, SR, and ST, and boundary terms from 
r = and r = oo given by 

J dt [\R 2 Q S^N^ 1 ) - M + SN + + Q 5% - Q + S$ + ) . (2.19) 

The variation thus gives the desired classical equations of motion provided we fix, in addition 
to the initial and final values of A, R, and T, also the quantities A^Aq , N + , $ , an d On 
a classical solution all these quantities have a clear geometrical interpretation. A'iAq 1 gives 
via ( |2.12j ) the evolution of the unit normal to the constant t hypersurface at the bifurcation 
two-sphere, and N + gives the evolution of the Killing time at the infinity. $o and $ + 
describe the electromagnetic gauge in a way that will become more transparent in Sec. [TV]. 
Note from ( |2.11i ) that when a classical solution is written in coordinates that are regular at 



the bifurcation two-sphere, the electromagnetic potential will be regular at the bifurcation 
two-sphere only if $ = 0. 



III. CANONICAL TRANSFORMATION 

In this section we perform a canonical transformation, which generalizes that given in 



Ref. [H for the spherically symmetric vacuum Einstein theory. Following Ref. we first 
examine how the variables appearing in the action (|2.17|) carry the information about the 
geometry of the classical solution ( |A1|) . We then use this information as a guide for finding 
the canonical transformation. 



A. Reconstruction 

Under our boundary conditions, every classical solution is an exterior region of a RNAdS 
spacetime with a nondegenerate Killing horizon (see Appendix |A|). We now assume that we 
are given the canonical data (A, R, T, P\, Pr, Pr) on a spacelike hypersurface embedded in 
such a RNAdS spacetime. We wish to recover from the canonical data the mass and charge 
parameters of the spacetime, the information about the embedding of the hypersurface in 
the spacetime, and the information about the electromagnetic gauge. 

Consider first the charge. The equations of motion imply that Pr is independent of both 
t and r. It is easily seen that in the curvature coordinates (|A1|) , the value of Pr is just the 
charge Q. As Pr is unchanged by the gauge transformations generated by the constraints, 
it follows that in any gauge 

Q = Pr • (3.1) 

Consider then the mass. The reconstruction of the function F appearing in the metric 
( |Ala| ) proceeds exactly as in Ref. 13], with the result 




Pa- 
R 



(3.2) 

From (| A 1 b| ) and ( |3.1| ), we find for the mass the expression 
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where F is understood to be given by (|3.2|). 

Consider then the embedding. By repeating the steps in Ref. fS^I , we obtain 



-T' = R^F^APa , (3.4) 

which determines the embedding up to an overall additive constant in T. To determine 
the value of the additive constant, one needs to know the value of T at one point on the 
hypersurface. 

Consider finally the electromagnetic gauge. By 
that 

A = R~ x QdT + 

= (R-^QT' + «£') dr + (R^QT + £) dt . (3.5) 
From (F3)> (|31]), and Q we then obtain 

£' = r + R- 2 F~ l AP A P T , (3.6) 
which determines the value of £ on the hypersurface up to an additive constant. 



A2), there exists a function £(t,r) such 



B. Transformation 

We have seen that when the equations of motion hold, the quantities defined by Eqs. 
( fj. 1| ) ( PT^D have a transparent geometrical meaning. We now promote these equations into 
definitions of functions on the phase space, valid even when the equations of motion do not 
hold. Our aim is to complete the set of functions into a set that constitutes a canonical 
chart. 

We shall from now on assume that the quantity R2 in Eq. (|2.11b|) is positive. As noted 
in Sec. [TT], this is always the case for our classical solutions. 



The functions M ( |3.3|) and Q ( |3 . 1|) Poisson commute with each other. The function —T' 
Poisson commutes with Q and is canonically conjugate to M. This suggests looking 
for a canonical transformation such that M and Q become two new coordinates, and — T" 
becomes the momentum conjugate to M. As in the Schwarzschild case |44||, the function 



R := R Poisson commutes with M, Q, and — T', and provides therefore a candidate for a new 
canonical coordinate. The crucial issue then is whether one can find momenta conjugate to 
Q and R such that the transformation is canonical. 

A necessary condition for the prospective new momenta arises from the observation that 
the supermomentum constraint ( |2.7b| ) generates spatial diffeomorphisms in all the variables. 



Since M, Q, and R are spatial scalars, the expression for H r in the new variables must be 
P M M' + P Q Q' + P R R'. Equating this with (p77B|) and substituting for M and P M = -V 
their expressions from ( |3.3|) and (|3.4j ) gives only one equation for the two unknowns Pq 



and Pr, but the structure of the equation as a linear combination of R' and Pp suggests 
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setting the coefficients of R' and P r individually to zero. These considerations suggest the 
transformation 



M 

Pm 
R 

Pr 



\R (r 



2/9-2 



+ 1 + P 2 R~ 2 - F 
'a , 



R^F^AP, 
R , 

P R - \R- l AP A - \R- l F~ 1 AP A 
-R^A^F- 1 ((AP A )'(PP') - (AP A )(RR 
+\R~ l F- 1 AP K [P^R~ 2 - 3P 



? 3«-2 



Q:=Pr 



P 



Q 



-r - R~ F~ AP\Pt 



(3.7a) 

(3.7b) 
(3.7c) 



(3.7d) 
(3.7e) 
(3.7f) 



where F is given by ( |3.2|) . The analogy between the pairs (M,P M ) and (Q,Pq) becomes 
manifest by observing from ( p.6| ) that on a classical solution, Pq carries the information 
about the electromagnetic gauge via Pq = — 

It is now straightforward to demonstrate that the transformation ( |3.7| ) is indeed canon- 
ical. We begin with the identity 



P A 5A+P R 8R + P r 5T - P M 5M - P R 5R - P Q 8Q 
RR' + AP A 



hR5R\n 



RR' - APa 



5 TP r + AP A + ±RR'\n 



RR' - AP A 



RR' + AP A 



(3.8) 



and integrate both sides with respect to r from r = to r = oo. The first term on the right 
hand side gives substitution terms from r = to r = oo that vanish by virtue of our falloff 
conditions, and we obtain 



/•OO 

/ dr (P A 5A + P R 5R + P r ST) 
Jo 

POO 

- / dr (P M 8M + P R 5R + P Q 5Q) = 5lu [A, P, T, P A , P r ] 
Jo " 



where 



uo [A, P, r, P A , P r ] = dr (vP r + AP A + \RR' In 



RR' - AP A 



RR' + APa 



(3.9) 



(3.10) 



The functional ui [A, P, T, Pa, Pr] is well-defined by virtue of the falloff conditions. Eqs. ( p.9| ) 
and ( p . 1 U| ) show that the Liouville forms of the old and new variables differ only by an exact 
form, and the transformation is thus canonical. 

The new variables have well-defined falloff properties at r = and r — > oo. At r = 0, 
Eqs. (|2~Tl| ) imply 



F(t,r) =4P^A"V + 0(r 



(3.11) 



and 
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where 



At r 



M(t, r) = M (t) + M 2 (t)r 2 + 0(r 4 ) , 


(3.12a) 


R(t,r) = Ro(t) + R 2 (t)r 2 + 0(r 4 ) , 


(3.12b) 


Q{t,r) = Q (t) + Q 2 (t)r 2 + O(r i ) , 


(3.12c) 


P M {t,r) = 0{r) , 


(3.12d) 


P R (t,r) = 0(r) , 


(3.12e) 


P Q (t,r) = 0(r) , 


(3.12f) 


M = \R, (R 2 r 2 + 1 + QlR- 2 ) , 


(3.13a) 


M 2 = \R 2 (3R 2 £- 2 + 1 - Q 2 Rf - 4RqR 2 A 2 ) + Q0Q2R0 1 ■ 


(3.13b) 


bqs. ( 2. ld|) imply 




M{t,r) = M+(t) + 0°°(r- 1 ) , 


(3.14a) 


R(t,r) = r + £ 2 p{t)r- 2 + 0°°{r- 3 ) , 


(3.14b) 


Q(t,r) = Q + (t) + 0(r- 1 ) , 


(3.14c) 


P M {t,r) = 0°°(r- 6 ) , 


(3.14d) 


P R (t,r) = 0°°(r- 4 ) , 


(3.14e) 


P Q (t,r) = 0°°(r- 2 ) , 


(3.14f) 



where M + (t) is given by ( ^.14| ). 

The canonical transformation (|3.7|) becomes singular when F — 0. Under our boundary 
conditions the classical solutions have F > for r > 0. At the limit r->0F approaches 
zero according to fl3.11|) , but ( |3.7|) has still a well-defined limit obeying (|3.12|) . Our canonical 
transformation is therefore well-defined and differentiable near the classical solutions, and 
similarly the inverse transformation is well-defined and differentiable near the classical so- 
lutions. From now on we shall assume that we are always in a neighborhood of the classical 
solutions such that F > holds for r > 0. 



C. Action 

It is possible to write an action in the new variables by simply re-expressing the con- 
straints (|2.7| ) in terms of the new coordinates and momenta. A more transparent action can 
be found if we exercise the freedom to redefine the Lagrange multipliers. 

The constraint terms in the bulk action (|2.6|) take the form 

NH + N r H r + $G = N M M' + N R P R + N Q Q' , (3.15) 

where 

N M = —NF~ 1 A~ 1 R' + N'R^F^APa , (3.16a) 

N R = -NR^Pa + N r R' , (3.16b) 

N Q = NR^F^A^R'Pr - N r (T + R^F^APaPA - I . (3.16c) 
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When viewed as a linear transformation from (JV, N r ,&) to (N M , N R , N®), Eqs. ( |3.16|) are 
nonsingular for r > 0. This suggests that we could take the constraint terms in the new bulk 
action to be those on the right hand side of ( |3.15| ), with N M , N R , and as independent 
Lagrange multipliers. At r — > oo this would be satisfactory: ( |3.16|) implies the asymptotic 
behavior 



N M (t,r) 
N R (t,r) 
N Q (t,r) 



-N + (t) + 
0°°(r- 2 ) 
-$+(*) + o 



oo/ r -l^ 



(3.17a) 
(3.17b) 
(3.17c) 



and one could then fix N+(t) and $+(t) as in Sec. [TT] after adding the boundary action 



- Jdt (N+M+ + $+Q. 



However, at r = we have 



N M (t,r) = -±N 1 A R^ 1 + O(r 2 ) , 
N R (t,r) = 0(r 2 ) , 

N Q (t,r) = -$ (t) + iNtAoQoR^ 1 !^ 1 + 0(r 



(3.18) 



(3.19a) 
(3.19b) 
(3.19c) 



which says that fixing N M and N® at r = to values that are independent of the canonical 
variables is not equivalent to fixing A^Ap 1 and $o to values that are independent of the 
canonical variables. We therefore need to redefine N M and N® near r = 0, without affecting 
their behavior at r — » oo. 

To proceed, we make two assumptions. First, we assume M > M crit (Q ), where the 
function M crit is defined in Appendix [A| Second, we regard Eq. (|3.13a|) as defining Rq in 
terms of M and Q as i? — -Rhor(A^o> Qo), where the function R hoT is defined in Appendix [A|. 
As discussed in Appendix [A], these assumptions are always true for our classical solutions, 
and they therefore merely tighten the neighborhood of the classical solutions in which the 
field variables may take values. For future use, we note that these assumptions imply 
2 > 0, and the variation of Rn takes the form 



3R 2 £- 2 + 1 - Q 2 Ro 



6Ro = 2{3R 2 £~ 2 + 1 - Q 2 Rf) 1 (SM - R l Q^Qv) 
Define now the quantities N M and A^ by 



N 



M 



-N 



M 



:i-g) + 2gR (3R 2 e- 2 + l-Q 2 R l 



-i 



(3.20) 

(3.21a) 
(3.21b) 

where g(r) is a smooth decreasing function that vanishes at r — >• oo as 0°°(r~ 5 ), and 
approaches the value 1 at r — > as g(r) = 1 + 0(r 2 ). Eqs. ( |3.21| ) then define a nonsingular 
linear transformation from (A^ M , N R ) to (A^ M , N®). The asymptotic behavior at r — > oo is 



n q = 2N M gQ {3R l 



QlR t 



2 d-2 




- N Q 



N M (t,r) 
N Q (t,r) 



N+(t) + 
$+(t) + 0°° 



(3.22a) 
(3.22b) 



14 



and the asymptotic behavior at r = is 



N M (t,r) =iV M (t) + O(r 2 ) 
N Q (t,r) = %(t)+0(r 2 ) . 



(3.23a) 
(3.23b) 



(3.24) 



where 

iV M = iJVxAoi^ 1 i? 2 - 1 (3i?gr 2 + 1 - QX~ 2X 

When the constraints M' = and Q' = hold, Eqs. fl3~T2aD , (gl2^) , and ftTIgED show that 

q'=o 



AT, 







(3.25) 



Thus, when the constraints hold, fixing N M and iV^ at r = is equivalent to fixing iViAg 1 
and $o- We therefore adopt N M , N R , and as a set of new independent Lagrange 
multipliers. 

The bulk action takes the form 



S*[M, R, Q, P M , P R , P q ; N m , iV R , iV Q ] 

= J dt J™ dr^P M M + P R R + P Q Q + N Q Q' - N R P R 



+N 



M 



(1 _ g) M ' + 2g{3R 2 £~ 2 + 1 - QlRf) ' (R M' - Q Q') 



(3.26) 



The total action is taken to be 



S[M, R, Q, P M , Pr, Pq; N m , N R , N Q ] = S S [M, R, Q, P M , P R , P Q ; iV M , N R , N Q ] 

+Sex[M ,M + ,Q ,Q + ;N +7 $ ,$ + ] , (3.27) 



where 



S as [M , M+, Q , Q+; N+, $o, $+] 



(3.28) 



The quantities to be varied independently are M, R, Q, Pm, Pr, Pq, N M , N R , and N®, 
and the boundary conditions for the new Lagrange multipliers are given by ( |3.17b| ), ( |3.19b ). 
( 13.22j ), and (|3.23|) . The volume term in the variation of ([3.271) is proportional to the equations 
of motion 



M = , 
R = iV R 
Q = , 



M 



N 



M 



(3.29a) 
(3.29b) 
(3.29c) 

(3.29d) 



15 



Pr = , 
Pq = (N Q )' 
M' = , 
Pr = , 
Q' = , 



(3.29e) 
(3.29f) 
(3.29g) 
(3.29h) 
(3.29i) 



where N M and are now understood to be defined by ( |3.21| ). The boundary terms 
in the variation consist of terms proportional to SM, 5R, and 6Q on the initial and final 
hypersurfaces, and terms from r = and r = oo given by 



dt 



M + 5N + + Q 5$o - Q+ 



(3.30) 



To arrive at (|3.30| ), ( |3.20| ) has been used. The action ( |3.27| ) thus yields the equations of 



motion ( 3.29|) provided that we fix, in addition to the initial and final values of the new 
canonical coordinates, also the quantities N^ 1 , N + , <3> , and $+. Because of ( [3.25D , these 
fixed quantities at the right and left ends have precisely the same interpretation in terms of 
the geometry of the classical solutions as the fixed quantities in the action (|2 . 1 7|) . 



IV. HAMILTONIAN REDUCTION 



In this section we shall reduce the action ( p.27|) to the true dynamical degrees of freedom 
by solving the constraints. 

The constraints (|3.29g|) and (|3.29i| ) imply that M and Q are independent of r. We can 
therefore write 



M(t,r) = m(t) 
Q(t,r) = q(t) 



(4.1a) 
(4.1b) 



Substituting this and the constraint (|3.29h|) back into ( |3.27|) yields the true Hamiltonian 
action 



S[m, q, p m , p q ; A^ M , A^ + , $ + , $ ] = / dt (p m m + p q q - h) 



where 



Pm 
Pq 



drP, 



M 



drPc 



The reduced Hamiltonian h in (1.2) is 

h = -\RlN^ + N + m + ($+ - $ ) q , 
where i?h := -Rhor(m, q). The assumptions made in the previous section imply 



(4.2) 

(4.3a) 
(4.3b) 

(4.4) 



16 



m > M crit (q) 



(4.5) 



and h is therefore well-defined. Note that h is, in general, explicitly time-dependent through 
the prescribed functions N^(t), N + (t), $+(£), and $ (*)- 

The variational principle associated with the reduced action (|4.2| ) fixes the initial and 
final values of the coordinates m and q. The equations of motion are 

m = , (4.6a) 
q = , (4.6b) 

p m = 2R h (3Rlr 2 + 1 - q 2 ^ 2 ) ^N^ 1 - N + , (4.6c) 
p q = -2q(3^ r 2 + l _ q 2^-2\ -iftM + $ _ $ + . ( 4 . 6d ) 



Eqs. ([4.6a ) and ( |4.6b| ) are readily understood in terms of the statement that on a classical 



solution m and q are respectively equal to the mass and charge parameters of the RNAdS 
solution. To understand Eq. (|4.6cj) , recall from Sec. [H| that on a classical solution P M = —T', 



where T is the Killing time. From ( [4.3a ) we see that p m = T — T + , where T and T + are 



respectively the values of T at the left and right ends of the constant t hypersurface. As the 
constant t hypersurface evolves in the RNAdS spacetime, the first and second term on the 
right hand side of ( |4.6c|) are respectively equal to To and — T+. The interpretation of Eq. 



Q4.6d|) is analogous. On a classical solution we have p q = £o — C+> where £ is the function 
that specifies the electromagnetic gauge via ( |3.5|) . The first two terms on the right hand 
side of (|4.6d|) give £o> an d the last term gives — £+. 

V. QUANTUM THEORY AND THE GRAND PARTITION FUNCTION 

We shall now quantize the reduced Hamiltonian theory of Sec. |V|. Our aim is to construct 
the time evolution operator in the Hamiltonian quantum theory, and then to obtain a grand 
partition function via an analytic continuation of this operator. 

A. Quantization 

As is well known, the quantization of a given classical Hamiltonian theory requires input 
fnH53|, and the questions of physically appropriate input for a quantum black hole remain 
largely open. For the purposes of the present paper we shall be content to define the quantum 
theory in essence by fiat, following Refs. Our main physical conclusions will emerge 



from the semiclassical regime of the theory, and at this level one may reasonably hope the 
details of the quantization not to be crucial. 

We regard m and q as configuration variables, satisfying the inequality ( |4.5|) . The wave 
functions are of the form ?/>(m, q), and the inner product is taken to be 

(^,X) = / udmdqipx , (5.1) 



A 
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where A C R 2 is the domain ( f4.5|) and /i(m, q) is a smooth positive weight factor. The 
Hilbert space is thus Ti := L 2 (A; [idmdq). We assume that \x is a slowly varying function, 
in a sense to be made more precise later, but otherwise it will remain arbitrary. 

The Hamiltonian operator h(t) is taken to act as pointwise multiplication by the function 
h(m, q; t) (|4.4|) : ip(m, q) i— > h(m, q; t)ip(m, q). h(t) is an unbounded essentially self-adjoint 
operator ||54j| , and the corresponding unitary time evolution operator in Ti is 



exp 



dt' h(t') 



k(t 2 ] ti) acts in Ti by pointwise multiplication by the function 



where 



K(m, q; T, S + , S , 6) = exp 



-imT 



zq(n+ 



H ) + URle 



t 2 



:-- 



*2 



t-2 



dtN+(t) , 
dt$+(t) , 
dt$ (t) , 
dtiV A/ (t) . 



(5.2) 

(5.3) 

(5.4a) 
(5.4b) 
(5.4c) 
(5.4d) 



k(t2',ti) therefore depends on t\ and t2 only through the quantities on the left hand 
side of (|5.4|) , and we may write K(t 2 ;ti) as K(T, H + , H , 0). The composition law, 
fT(t3; tzjKfa] t\) = K(ts; ti), amounts to independent addition in each of the four parame- 
ters in K(T, H + , Sq, Q), and we may regard these four parameters as independent evolution 
parameters specified by the boundary conditions. T is the Killing time elapsed at the in- 
finity, and 6 is the boost parameter elapsed at the bifurcation two-sphere. S + and S can 
be computed from the line integral of the electromagnetic potential ( |2.2| ) along the timelike 
curve of constant r and constant angular variables as this curve approaches respectively the 
infinity and the bifurcation two-sphere. 



B. Grand partition function 

We shall now construct a grand partition function by continuing the time evolution 
operator to imaginary time and taking the trace. We begin by discussing the boundary 
conditions for the relevant thermodynamical ensemble. 

The envisaged semiclassical thermodynamical situation consists of a charged spherically 
symmetric black hole in asymptotically anti-de Sitter space, in thermal equilibrium with a 
bath of Hawking radiation. If the back-reaction from the radiation is neglected, the geometry 
is described by the RNAdS metric ( |A1| ). Assuming that the local temperature is given in 



the usual manner in terms of the surface gravity and the redshift factor [|T8| , |3~5| , we see that 
the local temperature is F -1 / 2 /? -1 , where F is given by ( |Alb| ) and 
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p:= 47LRo(3r 2 f2jj + 1~Q 2 R 2 ) 1 . (5.5) 

At the infinity the local temperature vanishes as f3~ 1 £R~ 1 (1 + 0°°(£ 2 R~ 2 )), and (3~ l can thus 
be extracted from the asymptotic behavior as the coefficient of the leading order term £R~ X . 
We shall follow Refs. |TB| , P5]| and regard f3~ l as a renormalized temperature at infinity. 

The electromagnetic variable with thermodynamic interest for us is the electric potential 
difference between the horizon and the infinity, in the curvature coordinates ( |A1| ) and in 
an electromagnetic gauge that makes A invariant under the Killing time translations. We 
denote this quantity by <\>. From (|A2|) it is seen that on a classical solution <p = QRq 1 ■ 

We shall consider a thermodynamical ensemble in which the fixed quantities are /3 and 4>. 
This data can be interpreted as that for a grand canonical ensemble, with <fi being analogous 
to the chemical potential [ff|J|p6ll . Our aim is to obtain a grand partition function Z(/3, 0) by 
continuing the time evolution operator of the Lorentzian Hamiltonian theory to imaginary 
time and taking the trace. 

The continuation of T is straightforward: comparing the definition of (3 to the falloff of 
N in (HjJD and to the definition ( ggj ), we are led to set T = —i/3. For the continuation of 



we choose = —2iri, motivated by the regularity of the classical Euclidean solutions as 
in Ref. ||32|| . We mentioned at the end of Sec. || that the regularity of the electromagnetic 



potential at the bifurcation two-sphere of the Lorentzian solutions requires $o — 0; similarly, 
requiring regularity of the electromagnetic potential at the horizon of the classical Euclidean 
solutions now leads us to set So = 0. Finally, recall that S + gives the constant r line integral 
of the electromagnetic potential (|2.2j) at the infinity. Comparing this to the definition of 0, 
we set S + = — Tcf) = i/3<f). We are thus led to propose for the grand partition function the 
expression 

Z(0, (f)) = Tr [K(-i/3, if3cj), 0, -2m)] . (5.6) 

As it stands, the trace in ( |5.6|) is divergent, but one can argue as in Refs. f32| , |3"4|| that a 
suitable regularization and renormalization yields the result 

Zren((3, (ft) — N I ^dmd<\ exp —(3 (m — q</>) + nR 2 , (5.7) 



where we have substituted for K the explicit expression ( |5.3|) . The normalization factor J\f 
may depend on £, but we shall assume that it does not depend on (3 or (f>. 

Provided the weight factor \x is slowly varying compared with the exponential in (|5.7j) , it 
is easy to verify, using the definition of Rh given after Eq. (f4.4p, that the integral in ( |5.7[ ) is 
convergent. Equation (|5.7| ) thus yields a well-defined grand partition function. Comparing 
with ordinary PVT systems (f) is now indeed seen to be analogous to the chemical 
potential, and the quantities m and q are respectively analogous to the energy and the 
particle number. We shall examine the thermodynamical properties of this grand partition 
function in the next section. 



VI. THERMODYNAMICS IN THE GRAND CANONICAL ENSEMBLE 



It is useful to change the integration variables in 
(i?h,q)- From (|A6|) we obtain 



from the pair (m, q) to the pair 
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m = ±i2 h (Rlr 2 + 1 + q 2 i?h 2 ) , (6.1) 
and the grand partition function takes the form 

Z ven ((3,(j))=M [ fidR h dqexp(-Q , (6.2) 

J A' 

where 

h(R h , q) := \(3R h (Rlt 2 + 1 + q 2 i?h 2 ) ~ Pfa ~ *K • (6-3) 

One may view J* as an effective action or a reduced action [24H26I . The integration domain 
A 1 is given by the inequalities 

< R h , (6.4a) 
q 2 < Rl (l + 3Rlr 2 ) , (6.4b) 

and the weight factor \x is obtained from \i by including the Jacobian |<9(m, q)/d(Rh, q)|. 
Note that because of fl6.4b|) , 7* remains finite as R^ — > 0. 

As \i is assumed to be slowly varying, we can estimate Z ren (j3, 0) by the saddle point 
approximation to (|6.2|) . For this, we need to find the critical points of i* in the interior 
of A'. 

When 4> 2 < 1 - ^tt 2 £ 2 /3~ 2 , J* has no critical points. When 1 - ^n 2 f[3~ 2 < <fi 2 < 1, the 
two critical points of i* are at 



q = q± := (f)R± . (6.5b) 

The lower signs do not give a local extremum, but the upper signs give a local minimum. 
In the limiting case 1 — ^n 2 l 2 f3~ 2 = (ft 2 , the only critical point is (_R{j~,q + ), but it is not a 
local extremum. Finally, when 2 > 1, the only critical point is (i?J,q + ), and it is a local 
minimum. Whenever the critical points exist, the value of i* at these points can be written 

as 

ft( ^ ) = ^)'(i-^-(f>'^) . (6 „ 

1 h ^ 1 - 2 + 3(i£) V 2 

As grows without bound in the noncompact directions in A', the global minimum 
can be found by examining J* at the critical points and on the boundary of A'. When 
<p 2 > 1 — n 2 £ 2 j3~ 2 , the global minimum is at the critical point (R^, q + ), and h(R^ , q + ) is 
negative. When 2 < 1 — 7T 2 £ 2 /3~ 2 , the global minimum is at R^ = = q, where i* vanishes. 
In the limiting case <fi 2 = 1 — 7r 2 £ 2 j3~ 2 , vanishes at (i?{J",q + ) and at R^ = = q, and is 
positive everywhere else. 

We thus see that for <p 2 > 1 — ti 2 £ 2 /3~ 2 , Z ren can be approximated as 

Z ren (/3,0)^Pexp[-/*( J R+,q+)] , (6.7) 
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where P is a slowly varying prefactor. The approximation becomes presumably progressively 
better with increasing \I*(R£, q + )|. For 2 < 1 — 7r 2 £ 2 /? -2 , the dominant contribution to Z ien 
comes from the vicinity of = = q, and the behavior of Z Ten depends more sensitively 
on the weight factor jl. 

These results for Z ren are consistent with what one would expect just from the existence 
of (Lorentzian) black hole solutions under fixing <ft and the renormalized inverse Hawking 
temperature f3 (|5.5Q . It can be verified that such solutions exist precisely at the critical points 
of i*: the values of m and q at these critical points are just the mass and charge parameters 
of the black hole. Further, the value of I* at a critical point is simply the Euclidean action 
of the corresponding Euclideanized black hole solution. When a unique classical solution 
exists, it dominates the grand partition function; when two distinct classical solutions exist, 
the grand partition function is dominated either by the larger mass classical solution or by 
no classical solutions. The situation is thus remarkably similar to that found in the absence 
of a cosmological constant when the boundary conditions are set on a finite size box ||26|| . 

Let us now consider the thermodynamical predictions from Z ren . Recall that the thermal 
expectation values of the energy and charge in the grand canonical ensemble are given by 

(E) = [~ + ^^) 0*2**) , (6.8a) 
(QH/r 1 ^^ . (6.8b) 

When Z ren is dominated by the critical point (R^, q + ), we find 

(E) « m+ , (6.9a) 
(Q)^q+ , (6.9b) 



where m + is obtained from (R^ , q + ) through (|6.1| ). That is, the thermal expectation values 
of the energy and the charge are simply the mass and charge parameters of the dominant 
classical solution. In particular, there are no additional contributions to the mass from the 
gravitational binding energy associated with the thermal energy, or from the electrostatic 
binding energy associated with the charge. Such additional, finite size contributions were 
found to be present in the finite size ensembles of Refs. [ II , 21 , 25 , 51 1. 



It is easily seen that (dm + /dj3) < 0. This means that when the approximation 
( |6.9a| ) is good, the (constant (f>) heat capacity, = —j3 2 (d(E)/dj3), is positive. In the 



regime ( |6.9a[ ), the system is thus stable under thermal fluctuations in the energy. Note 
that as (dm~/d/3) > 0, a grand partition function dominated by the lower mass classical 



solution would be thermodynamically unstable This is analogous to what happens in 
the absence of a cosmological constant under the boxed boundary conditions considered in 

Refs. ran 



It is also easily seen that (<9q + /<90) > 0. This shows that when the approximation ( |6.9b|) 
is good, we have (d(Q)/d<j>) > 0, and the system is stable under thermal fluctuations in the 
charge. More generally, one can show directly from the expressions ( |6.2j ), ( |6.3| ), and ( |6.8b| ) 
that (d(Q)/d(p) > holds always, even when the approximation (|6.9b|) is not good. 

The entropy in the grand canonical ensemble is given by 



21 




S = U-fe (In^ren) • (6-10) 



When the approximation (]6.7|) is good, we have 5 ~ 7r(i?^) 2 , which means that the entropy 



is one quarter of the horizon area. This is the anticipated Bekenstein-Hawking result. 

Finally, when Z Ten is not dominated by a critical point, the thermodynamical predictions 
become much more sensitive to the choice of the weight factor jx. As in Refs. p^- p6| . |34| , one 
can view the transition in the qualitative behavior of Z ien as evidence for a phase transition 
between a black hole sector and a topologically different sector of the theory; in the case 
at hand, the second sector might be referred to as "hot anti-de Sitter space." On classical 
grounds one might have expected this transition to occur near <fr 2 « 1 — ^7i 2 £ 2 /3~ 2 , where 
the classical solutions disappear. However, we saw that the transition in fact occurs near 
<p 2 « 1 — 7T 2 £ 2 (3~~ 2 , where the two classical solutions still exist. This is highly similar to 
what happens in four dimensions under boxed boundary conditions without a cosmological 
constant |2"1| , |2"6] | , but subtly different from what happens in two dimensions with Witten's 
dilatonic black hole [M . 



VII. THE CANONICAL ENSEMBLE 

We have seen that the Hamiltonian formulation of Sees. [I I - [TV| led into a thermody- 
namical grand canonical ensemble where the fixed quantities are the renormalized inverse 
temperature (3 at infinity and the electric potential difference between the horizon and the 
infinity. From the thermodynamical viewpoint, another natural ensemble for the charged 
black hole in asymptotically anti-de Sitter space is the canonical ensemble, where one allows 
fluctuations in but fixes instead the charge q. In this section we shall outline the recovery 
of the canonical ensemble from a Lorentzian Hamiltonian analysis, and briefly discuss the 
thermodynamical properties of the black hole in this ensemble. 

As a starting point, we modify the boundary conditions of the Hamiltonian theory 
of Sec. [TJ by leaving $o(t) and unspecified but fixing Qo(t) and Q+(t) to be pre- 

scribed functions of t. The action is obtained from (|2.17 ) and (|2.18|) by omitting the terms 



/ dt (&0Q0 — from fl2.18|) . Clearly, classical solutions exist only when Qo(t) and 



Q+(t) are chosen independent of t and equal. We shall from now on assume that the bound- 
ary data is chosen in this manner. 

One way to proceed is simply to push through the canonical transformation of Sec. 
noting that the new boundary conditions merely result into minor modifications. It is only 
when one subsequently performs a Hamiltonian reduction along the lines of Sec. [IV] that the 
new boundary conditions give rise to important differences. Firstly, the boundary data for 
Qo(t) and Q+(t) implies that the quantity q(i) defined by ( |4.1b| ) is a t- independent constant 
whose value is completely determined by the boundary conditions. Therefore, the Liouville 
term / dt p q q drops entirely out of the action (|4.2j ). Secondly, because of the terms that were 
omitted from the boundary action ( |2.1£| ), the term (& + — $ ) Q drops out of the reduced 
Hamiltonian ( (4.4|) . This means that in the reduced Hamiltonian theory q has become an 
external parameter specified by the boundary conditions: it is not varied in the action, and 
it does not have a conjugate momentum. The new reduced action reads 
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Sb[m, p m ; N M , N + ; q] = J dt (p m m - h c ) , (7.1) 

where 

h c = -±RlN^ + N+m . (7.2) 

Here Ph := Phor( m , q) as before, and the assumptions made in the canonical transformation 
again imply that ( f4.5|) holds. 

An alternative way to proceed under the new boundary data is to partially reduce the 
action already in the variables of Sec. || by solving the constraint (|2.9c|) . One uses the 
constraint ( |2.9c| ) and the equation of motion ( |2.101| ) to set Pr(t, r) equal to the constant 



specified in the boundary data, and substitutes this back in the action. The Liouville term 
J °° dr PpT then becomes a total time derivative and can be dropped. One thus obtains an 
action that no longer involves T or $, involves Pp only as a prescribed constant, and correctly 
yields the equations of motion for the remaining variables. One can now perform a canonical 
transformation from the variables (A, R, Pa, Pr) to the new variables (M, R, P M , P R ), defined 



as in Sec. [I I J except that Pr = Q is now regarded as a fixed external parameter. Finally, 
one can reduce the action by solving the constraints as in Sec. |V[ The result is again the 
action given by ( |7.1|) and (|7T^ ). 

Quantization of the reduced Hamiltonian theory proceeds as in Sec. 0. For the renor- 
malized trace of the analytically continued time evolution operator, we obtain 



Z ren (/5,q)= J fidR h exp (-I c *) , (7.3) 

■Rcrit(q) 

where the function P cri t is defined by Eq. ( |A4|) in Appendix the weight factor jl is a 
positive function of Ph (and possibly q), and 

Ic*(Rh) ■= iPRh (Rlr 2 + 1 + q 2 Ph 2 ) - vrP 2 . (7.4) 

Under the assumption that jl is slowly varying, the dominant contribution to Z rcn can 
be estimated by saddle point methods. The cases q = and q ^ merit each a separate 
analysis. 



Consider first the special case q = 0. The lower limit of the integral in ( |7.3| ) is then at 
Ph = 0. The critical point structure of Ic* is identical to that of I* ( [7.4|) for (p = 0, and the 
locations of the critical points and the values of the action at these points can simply be 
read off from Sec. by setting = 0. 

Consider from now on the generic case q ^ 0. Ic* has one negative critical point, and 
from one to three positive critical points. The negative critical point is unphysical, but all 
the positive critical points lie in the physical domain P h > P crit (q). As Ic* is decreasing at 
Ph = Pcrit (q) and tends to infinity as Ph — > oo, the global minimum of Ic* in the domain 
Ph > P cr i t (q) is at a critical point. We can therefore concentrate on the positive critical 
points. 

When f3 2 > ^n 2 £ 2 , Ic* has only one positive critical point. When f3 2 < ^ir 2 £ 2 , the number 
of positive critical points is determined by the status of the double inequality 
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1 -"')(! + ') < ^ < HiMM , (7.5) 



36(1 -s) z ~ ~ 36(1 + s 

where 



2/32 (7.) 



3tt 2 £ 2 



When ( |7. 5| ) does not hold, there is only one positive critical point. When ( [7.5| ) holds as a 
genuine inequality, there are three positive critical points, and saturating the inequalities 
gives limiting cases where two of the three positive critical points merge. (Note that if 
(3 2 < f tt 2 ^ 2 , the leftmost expression in (|7.5|) is non-positive, and the left hand side inequality 
is then necessarily genuinely satisfied.) Now, when only one positive critical point exists, 
this critical point is the global minimum. On the other hand, when three positive critical 
points exist, they constitute a local maximum between two local minima, and the global 
minimum can be at either of the local minima depending on the values of the parameters. 
For example, when the right hand side inequality in (|7.5|) is close to being saturated, the 
global minimum is at the local minimum with the larger value of R^. 

The critical points can be examined further by parametrizing (3 and q as 



= („-») ,[3 (u W) + 1| 
u -\- v — uv 



(7.7a) 



q2r2 = «vp» , 

u -\- v — uv 

where the parameters u and v satisfy < v < u. The negative, unphysical critical point is 
then at = —£v, and R^ = iu gives a positive critical point.Q The condition that only one 
positive critical point exist reads 

a(u) < v , (7.8) 

where a(u) is the unique solution to the equation 

= 9ua 3 - (6m 2 + l) a 2 + u (9m 2 + 2) a - u 2 (7.9) 

in the interval < a < u. In this case the parametrization ( [7.7| ) is unique. When the 
inequality in (|7.8|) is reversed and three positive critical points exist, the parametrization 
( [7.7[ ) can be made unique by imposing the conditions 

«<-/! > ( 7 - 10a ) 



v ;1 + 6m2)(1_3m2)-(1-3m 2 ) 
v < 1 , (7.10b) 



*We thank Bernard Whiting for suggesting this type of parametrization. 
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which make i?h = £u the local maximum. The two local minima are then at the roots of the 
quadratic equation 

= 3(V + v 2 -uv)(R h /£) 2 - {u-v){3uv + l)(R h /£) +uv(3uv + l) . (7.11) 
The global minimum is at the larger (smaller) local minimum when the inequality 

< 12(6w- l)(u 2 + v 2 -uvf + (u-v) 2 (3uv + 1) [3 (u 2 + v 2 ) + l] (7.12) 
is satisfied (reversed). 

It is of some interest to examine the behavior of the critical points in the limit q 2 — > 
with fixed j3. When j3 2 > ^ti 2 £ 2 , the above discussion shows that for sufficiently small q 2 
there exists only one positive critical point, and in the limit q 2 — > this critical point 
approaches zero as 

R h = |q| [l + 2vr/r 1 |q| + O (q 2 *~ 2 )l . (7.13) 



When P 2 < |7r 2 £ 2 , on the other hand, there exist three positive critical points for sufficiently 
small q 2 . In the limit q 2 — > 0, the smallest positive critical point again approaches zero 



as ( [7.13|) , whereas the two larger ones approach the two critical points of the case q = 0. 



In the limiting case /3 2 = |7r 2 £ 2 , the smallest of the three positive critical points once again 
approaches zero as ( |7.13| ), and the two larger ones merge into a q = critical point that 
is not a local extremum. The limiting behavior is thus smooth, in spite of the changing 
number of critical points. 

At any critical point, the value of the action can be written as 

ttR 2 to' 2 - 1 - 3q 2 i? h 2 ) 
c * 3R 2 e~i + l-cfR^ ' 1 } 

In the limit q — > 0, this agrees with the expression given in Ref. [fJIJ. 

We thus see that for generic values of the parameters, Z Ten (fi, q) can be approximated 
by exp —Icl n > where /™™ stands for the value of Ic* at the critical point that is the global 
minimum. As in Sec. [Vl|, this is consistent with what one would have expected just from the 
existence of (Lorentzian) black hole solutions under fixing the charge and the renormalized 
inverse Hawking temperature: such solutions exist precisely at the critical points of Ic*, 
and the values of m and q at these critical points are just the mass and charge parameters 
of the black hole. One may view the shifting of the global minimum of Ic* from one local 
minimum to the other as a thermodynamical phase transition. 

We end this section with some brief remarks on the thermodynamics in the canonical 
ensemble. Recall that the formulas for the thermal expectation values for the energy and 
the electric potential read 

(E) = ^— , (7.15a) 

W = -/r 1 ^^ ■ (r.isb) 
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When a critical point of Ic* dominates, we obtain 



(E) « m , (7.16a) 
. (7.16b) 

These are, respectively, just the mass and the electric potential difference between the hori- 
zon and the infinity for the dominating classical solution. When the approximation (|7.16a ) 



is good, the positivity of the (constant q) heat capacity, C q = — (3 2 (d(E) /dj3), follows from 
the fact that the dominant critical point is a minimum of Ic* |24|j2~5|| . The positivity of C q 
follows more generally, even when the saddle point approximation does not hold, by direct 
manipulations from the expression ( |7.3|) and the assumption that ft is positive. 

When the saddle point approximation is good, we have for the entropy the Bekenstein- 
Hawking result, S = (1 - f3(d/d{3)) (lnZ ren ) « n(R h ) 2 . 

VIII. CONCLUSIONS AND DISCUSSION 

In this paper we have investigated the Hamiltonian dynamics and thermodynamics of 
spherically symmetric Einstein-Maxwell theory with a negative cosmological constant. We 
first set up a classical Lorentzian Hamiltonian theory in which the right end of the space- 
like hypersurfaces is at the asymptotically anti-de Sitter infinity in an exterior region of a 
RNAdS black hole spacetime, and the left end of the hypersurfaces is at the bifurcation two- 
sphere of a nondegenerate Killing horizon. We then simplified the constraints by a canonical 
transformation, and we explicitly reduced the theory into an unconstrained Hamiltonian the- 
ory with two canonical pairs of degrees of freedom. The reduced theory was quantized by 
Hamiltonian methods, and a grand partition function for a thermodynamical grand canoni- 
cal ensemble was obtained by analytically continuing the Schrodinger picture time evolution 
operator to imaginary time and taking the trace. The analytic continuation at the bifur- 
cation two-sphere was done in a way motivated by the smoothness of Euclidean black hole 
geometries as in Ref. [Q. A similar analysis with minor modifications to the boundary 



conditions led to a partition function for a thermodynamical canonical ensemble. Both the 
canonical ensemble and the grand canonical ensemble turned out to be well defined, and we 
were able to find the conditions under which the (grand) partition function is dominated by 
a classical Euclidean solution. 

Both thermodynamical ensembles exhibited a phase transition. In the grand canonical 
ensemble the transition occurs when the grand partition function ceases to be dominated 
by any classical Euclidean black hole solution, in close analogy with what happens in the 
spherically symmetric vacuum canonical ensemble with a finite boundary |I^ , p4| , p5| . In 
the canonical ensemble this kind of a phase transition can occur only in the limit of a 
vanishing charge, whereas for nonvanishing charge there occurs a phase transition in which 
the dominating contribution to the partition function shifts from one classical Euclidean 
solution to another as the boundary data changes. In either ensemble, whenever the (grand) 
partition function is dominated by a classical solution, one recovers for the entropy the 
Bekenstein-Hawking value of one quarter of the horizon area. 
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The classical canonical transformation of Sec. [TTT| is a relatively straightforward general- 
ization of the transformation that was found by Kuchaf in the spherically symmetric vacuum 
Einstein theory under Kruskal-like boundary conditions JHj]. When the classical equations 
of motion hold, our new canonical coordinates M and Q are simply the mass and charge 
parameters of the RNAdS solution. By (generalized) Birkhoff's theorem, the spacetime is 
uniquely characterized by these two parameters and the cosmo logical constant. The con- 
jugate momenta, Pm and Pq, carry the information about the embedding of the spacelike 
hypersurface in the spacetime and the electromagnetic gauge. Upon elimination of the con- 
straints, we saw in Sec. [TV] that Pm and Pq each give rise to one unconstrained momentum in 
the reduced Hamiltonian theory. These reduced momenta are global constructs with no local 
geometrical meaning, and they are associated with the anchoring of the spacelike hypersur- 
faces at the infinity and at the bifurcation two-sphere. The electromagnetic pair (Q, Pq) is 
quite closely analogous to the gravitational pair (M, Pm)- The third canonical pair, (R, P R ), 
is entirely gauge, and it completely disappears when the constraints are eliminated. 

Although we have here discussed the canonical transformation only under boundary con- 
ditions motivated by our thermodynamical goal, it would appear possible to use arguments 
similar to those in Refs. [0,0] to adapt this canonical transformation to boundary condi- 
tions under which the spacelike hypersurfaces extend from a left hand side asymptotically 
anti-de Sitter region to a right hand side asymptotically anti-de Sitter region, crossing the 
event horizons in arbitrary ways. The form ( 3.15Q taken by the constraints then suggests 
that, after introducing electromagnetic variables analogous to the reparametrization clocks 
t± of Ref. ||44j| , it is possible to perform a canonical transformation that separates Q into 



the charge density Q' and the charge at the (say) left hand side infinity, in analogy with the 
transformation that in Ref. [44] separates M into the mass density M' and the mass at the 
left hand side infinity. Also, it appears possible to take the limit where the cosmological con- 
stant vanishes and the asymptotically anti-de Sitter regions are replaced by asymptotically 
flat regions.^ It would further be possible to consider boundary conditions of the kind put 
forward in Ref. ||57|| . We have not investigated these issues in a systematic fashion; however, 
we shall outline in Appendix |B| how our canonical transformation can be adapted to the 
limit of a vanishing cosmological constant, under boundary conditions that still keep the left 
end of the hypersurfaces at the bifurcation two-sphere of a nondegenerate Killing horizon 
but replace the asymptotically anti-de Sitter falloff conditions at the right end by asymp- 
totically flat falloff conditions. In this case, each classical solution is the exterior region of a 
non-extremal Reissner-Nordstrom black hole. 

The thermodynamical results of Sees. [V4| and [VII] show that the stabilizing effect of the 
negative cosmological constant is highly similar to the stabilizing effect of a finite "box" with 
fixed surface area and fixed local temperature |l4] . |2~4] - f^ , |34]l . One important difference is, 
however, that in the asymptotically anti-de Sitter case various thermal expectation values 



5 Before the work reported in this paper was begun, we were informed by Karel Kuchaf that he 
had generalized the canonical transformation of Ref. [44] to the spherically symmetric Einstein- 
Maxwell system without a cosmological constant fl56| . We thank Karel Kuchaf for correspondence 
on this point. 
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are more directly related to the parameters of the dominant classical solutions. In the grand 
canonical ensemble, equations (|6.9| ) show that the thermal expectation values of energy and 
charge are simply the mass and charge parameters of the dominant classical solution: there 
are no additional contributions to the mass from the gravitational binding energy associated 
with the thermal energy, or from the electrostatic binding energy associated with the charge. 
Such additional, finite size contributions were found to be present in the finite size ensembles 
considered in Refs. [0,^,^,^]. In the canonical ensemble, the situation is similar with the 
thermal expectation values of the energy and the electric potential ( |7.16| ). 

The stabilizing effect of the negative cosmological constant becomes fully apparent when 
one attempts to repeat the analysis with a vanishing cosmological constant, replacing the 
asymptotically anti-de Sitter infinity by an asymptotically flat infinity. We shall outline this 
analysis in Appendix |B|. While there is no difficulty in quantizing the reduced Hamiltonian 
theory, the trace of the analytically continued time evolution operator turns out to remain 
divergent even after a renormalization of the kind performed in Sees. V B and |V11| . Neither 
the canonical ensemble nor the grand canonical ensemble exists. For the canonical ensemble 
this conclusion might be surprising in view of the observation that a Reissner-Nordstrom 
black hole in asymptotically flat space is stable against Hawking evaporation when one fixes 
the charge and the temperature at the infinity, provided the mass and charge parameters 
of the hole satisfy the inequality q 2 > |m 2 0. However, as we shall see in Appendix |B], 
the local stability of a classical solution is not sufficient to guarantee the existence of a full 
thermodynamical ensemble. 

Finally, we recall that as the physical temperature of Hawking radiation is redshifted 
to zero at the anti-de Sitter infinity, we followed Refs. and defined a renormal- 

ized temperature at infinity in terms of the rate at which the local Hawking temperature 
approaches zero. This definition led to physically reasonable conclusions; in particular, we 
recovered from the thermodynamical ensembles the Bekenstein-Hawking result for the black 
hole entropy. The definition can however be argued to have an ad hoc flavor, and one might 
wish to replace it by something that can be given a more immediate physical justification.^ 
What would be needed is a better understanding as to whether asymptotically anti-de Sitter 
infinity can in some sense be regarded as a physically realizable system, rather than just as 
a mathematically elegant set of boundary conditions. 

Note added. After this work was completed, we became aware of Refs. |58|-|60|j, which 
discuss the Dirac quantization of four-dimensional spherically symmetric Einstein-Maxwell 
geometries and related dilatonic theories. The work in these references has close technical 
similarities to our work. 



6 An interesting possibility might be to assume the Bekenstein-Hawking entropy and then to derive 
the appropriate renormalized temperature [38]. However, it does not appear clear how to adopt 
this as a starting point in a theory where the Bekenstein-Hawking result is expected to be only an 
approximate one, in the domain where the (grand) partition function is dominated by a classical 
Euclidean solution. 
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APPENDIX A: REISSNER-NORDSTROM- ANTI-DE SITTER BLACK HOLE 



In this appendix we recall some relevant properties of the Reissner-Nordstrom-anti-de Sit- 
ter (RNAdS) metric. We concentrate on the case where a nondegenerate event horizon exists, 
and on the region exterior to this horizon. 

In the curvature coordinates (T, R), the RNAdS metric is given by 

ds 2 = —FdT 2 + F~ l dR 2 + R 2 dQ? , (Ala) 
where dfl 2 is the metric on the unit two-sphere and 

R 2 2M Q 2 , . 

F -= T + 1 -1T + W ' < Alb > 

T and R are called respectively the Killing time and the curvature radius. The parameter £ is 
positive, and we take the parameters M and Q to be real. Together with the electromagnetic 
potential 

A = % dT » (A2) 



the metric (|A1|) is a solution to the Einstein-Maxwell equations with the cosmological con- 



stant — 3£~ 2 pE0| , f43f . The parameters M and Q are referred to respectively as the mass and 
the (electric) charge. The case Q = yields the Schwarzschild-anti-de Sitter metric, and the 
case Q = M = yields the metric on (the universal covering space of) anti-de Sitter space 

The metric ([AT]) has an asymptotically anti-de Sitter infinity at R — > oo for all values 



of the parameters [47j. We wish to restrict the parameters so that the metric describes 
the exterior of a black hole with a nondegenerate horizon. This happens when the quartic 
polynomial R 2 F(R) has a simple positive root R = Ro, such that F is positive for R > Rq. 
The necessary and sufficient condition is M > M CTit (Q), where 

M crit (Q) := ^= (yi + 12(Q/£) 2 + 2) (y/l + \2{Q/lf - l) ^ . (A3) 

Note that M is then necessarily positive. Rq can now be determined uniquely as the function 
Rhor{M, Q) of M and Q: for Q = 0, Rh or (M, Q) is defined as the unique positive solution 
to the equation F = 0; for Q 7^ 0, Rh OT (M, Q) is defined as the larger of the two positive 
solutions. In either case, if Q is considered fixed, Rhor(M, Q) is a monotonically increasing 
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function of M that takes the values R cr it(Q) < Rhor(M,Q) < oo as M cr i t (Q) < M < oo, 
where 



R CTit (Q) : = — [^l + 12(Q/£f-l 



1/2 



(A4) 



The metric can thus be uniquely parametrized by Q and Ro- The only restriction for 
these parameters is 



-Ro > -Rcrit(Q) ? 



(A5) 



and the mass is then given by 

^t(I + 1 + f) • < A6 > 

With Rq < R < oo, the metric ( |Ala| ) covers the region from the horizon to the asymp- 



totically anti-de Sitter infinity. The Penrose diagram can be found in Refs. [37,52 



APPENDIX B: REISSNER-NORDSTROM BLACK HOLE IN 
ASYMPTOTICALLY FLAT SPACE 

In the main text we took the cosmological constant to be strictly negative. In this 
appendix we shall outline the corresponding classical and quantum mechanical analysis in 
the case where the cosmological constant vanishes. In the notation of the main text this 
means taking the limit £ — ■> oo. The classical solutions are then not asymptotically anti- 
de Sitter but asymptotically flat, and the falloff conditions at r — ► oo must be modified to 
reflect this fact. 

In the variables of Section [TT], we retain the falloff conditions (|2.11 ) at r —>■ 0, but at 
r — > oo we introduce the new falloff conditions 



A(t,r) 


= l + M + (t)r- 1 + O co (r- 1 - e ) , 


(Bla) 


R(t,r) 


= r + 0°°(r~ e ) , 


(Bib) 


PA(t,r) 


= 0°°(r- e ) , 


(Blc) 


P R (t,r) 


= O co (r- 1 - e ) , 


(Bid) 


N(t,r) 


= iV + (t) + 0° o (r-) , 


(Ble) 


N r (t,r) 


= 0°°(r- e ) , 


(Blf) 


T{t,r) 


= O co (r- 1 - e ) , 


(Big) 


P r {t,r) 


= Q + (t) + 0°°(r- e ) , 


(Blh) 


r) 


= $ + (t) + 0° o (r-) , 


(Bli) 



where < e < 1. For the metric quantities these conditions are precisely those used in Ref. 
44| , ensuring asymptotic flatness. These conditions make the bulk action (2^) well-defined, 



and they are preserved under the time evolution. Adding the boundary action 
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J dt (|^iViA 1 - N + M + + $ Qo - (B2) 

yields an action for a variational principle in which N + , NiAq 1 , $ + , and $o are prescribed 
functions of t. Dropping the last two terms in ( p2|) yields an action for a variational principle 
in which $ + and $o are free but Q + and Qq are prescribed. 

The canonical transformation of the main text can now be adapted to the present bound- 
ary conditions by simply taking the limit £ — ► oo. A new action can be constructed as in 



Sec. [Ill Q the new falloff conditions only giving rise to minor technical modifications to the 
redefinition of the Lagrange multipliers. In the theory that prescribes $ + and $o ; elimination 
of the constraints along the lines of Sec. fV| yields the reduced action 

3[m, q, p m , p q ; iV M , N + , $+, $„] = J dt (p m m + p q q - h) , (B3) 

where the reduced Hamiltonian is given by 

h = -\RlN™ + iV+m + ($+ - $ ) q (B4) 



with R^ := m + y/m 2 — q 2 . The range of the variables is < m, q 2 < m 2 . In the theory 
that prescribes Q + and Q , one proceeds as in Sec. |V11| to obtain the reduced action 

S c [m, p m ; iV M , N + ; q\ = J dt (p m m - h c ) , (B5) 

where q is now regarded as an external parameter and 

h c = -lR 2 h N™ + N + m . (B6) 

Quantization of the two reduced theories proceeds as in the main text. For the renor- 
malized trace of the analytically continued time evolution operator, we obtain formally 

Z mn {l3,(t))=U J j2dR h dqexp(-Q , (B7a) 
Rh>h\ 

oo 

Zren(f3, q) = / fidR h exp (-I c *) , (B7b) 
|q| 

where I* and Ic* are respectively given by dropping the term proportional to £~ 2 from 
Eqs. (|6.3j) and (|7.4|). (3 is now interpreted as the inverse Hawking temperature at the 



infinity, with no renormalization. However, both integrals in (p7| ) are divergent because 
of the behavior of I* and Ic* at large R^. Thus, neither the canonical ensemble nor the 
grand canonical ensemble exists under the asymptotically flat boundary conditions. In this 
respect, the inclusion of the charge has therefore not made a qualitative difference from the 
asymptotically flat vacuum case p2 |. 



The critical points of I* and Ic* give again the (Lorentzian) classical solutions that have 
the inverse Hawking temperature (3 at infinity and the prescribed value of respectively <fi 
or q. The condition that I* possess critical points is \<p\ < 1: when this condition is satisfied 
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there exists exactly one critical point, but this critical point is not a local extremum. This 
reproduces the observations made by Davies in Ref. about charged black hole equilibria 
with fixed 0, and reflects in particular the fact that a semiclassical charged black hole under 
these boundary conditions is not stable against Hawking evaporation. 

The condition that Ic* possess critical points is /3/|q| > 6tc\/3, and when the inequality 
is genuine, there exist two critical points. The critical point with the smaller (larger) value of 
i?h is a local minimum (maximum, respectively). The local minimum satisfies q 2 > f m 2 , and 
it corresponds to the classical solution that Davies [|7| showed to be stable against Hawking 
evaporation under these boundary conditions (see also Refs. While the thermody- 

namical stability of this semiclassical solution is reflected in its being a local minimum of 
our Ic* |25| , the divergence of the the integral in ( |B7b| ) demonstrates that this local stability 
is not sufficient to guarantee the existence of a full thermodynamical canonical ensemble. 
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